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We consider an ideal gas of massive hadrons in thermal and chemical equilibrium. The finite-size 
gas expands longitudinally in accordance with Bjorken law. Also the transverse expansion in a form 
of the rarefaction wave is considered. We show that J/ty suppression in such an environment, when 
combined with the disintegration in nuclear matter, agrees qualitatively well with NA38 and NA50 
data. 



I. INTRODUCTION 

Since the paper of Matsui and Satz 0j there is a steady interest in the problem of J /\P suppression in a heavy-ion 
collision. The question is if this suppression can be treated as a signature for a quark-gluon plasma or if it can be 
explained by J/'F absorption in a hadron gas which appears in the central rapidity region (CRR) of the collision [§-|. 

In the following paper we shall continue our previous investigations |6| into the problem of J suppression observed 
in a heavy- ion collision (for experimental data see e.g. j?| and references therein). Now, we shall focus on the 
dependence of the suppression on the initial energy density reached in the CRR. 

In our model, J/^P suppression is the result of a cc state absorption in a dense hadronic matter through interactions 
of the type 

cc + h — >D + D + X, (1) 

where h denotes a hadron, D is a charm meson and X means a particle which is necessary to conserve the charge, 
baryon number or strangeness. The hadronic matter is in a state of an ideal gas of massive hadrons in thermal and 
chemical equilibrium and consists of all species up to baryon. Time evolution is given here by conservation laws 
combined with assumptions about the space-time structure of the system. A corresponding equation of state of the 
ideal gas makes then possible to express gas parameters such as temperature and chemical potentials as functions of 
time. 

An ideal gas of real hadrons has a very interesting feature: it cools much slower than a pion gas when expands 
longitudinally. We have checked numerically that for the initial energy densities eo corresponding to initial temper- 
atures To of the order of 200 MeV and for the freeze-out 7/. . not lower than about 100 MeV, the time dependence 
of the temperature of the expanding gas still keeps the well-known form T(t) — To ■ t~ a (we put to = Only 
the exponent a changes from ^ for massless pions to the values 5^-5^3 for massive realistic hadrons. As a result, 
the time of the freeze-out i/. G is much greater for the hadron gas than for the pion one. For instance, when we 
take T = 200MeV and T f , Q . = UOMeV we obtain t f . . = 7.37/m (a = for the hadron gas and t f . . = 2.9 fm 
(a = |) for the pion gas. The lower 1/ , the stronger difference. For T = 200MeV and Tf = lOOMeV we have 
i/. . = 48.5/m (a = ^g) for the hadron gas and tf. . = 8fm (a = ^) for the pion gas. This has a direct consequence 
for J suppression: the longer the system lasts, the deeper suppression causes (see Fig. |l|). 
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FIG. 1. comparison of suppression of the pure J/^'s in the cooling hadron gas for two values of the power a in approximation 
T(t) To • t~ a : a) a = | ; b) a = -L 

We are going to calculate a survival factor for J/ty when new, more realistic conditions are taken into account. 
We consider a hadronic gas which is produced in the CRR region. This gas expands both longitudinally and trans- 
versely. The longitudinal expansion is a traditional adiabatic hydrodynamical evolution Q , the transverse expansion 
is considered as the rarefaction wave. An initial energy density depends now on the impact parameter b and on the 
geometry of the collision. 



II. THE EXPANDING HADRON GAS 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of I species of particles, energy density 
e, baryon number density ns, strangeness density ns and entropy density s read {% = c = 1 always) 

— > (2 Si + 1) / r — , (2a) 
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where = (m 2 + p 2 ) 1 / 2 and m^, Sj, /ij, s, and <?j are the mass, baryon number, strangeness, chemical 
potential, spin and a statistical factor of specie i respectively (we treat an antiparticle as a different specie). And 
Hi = BtfiB + Si [is, where n B and n s ar e overall baryon number and strangeness chemical potentials respectively. 

We shall work here within the usual timetable of the events in the CRR of a given ion collision (for more details 
see e.g. ||). We fix t = at the moment of the maximal overlap of the nuclei. After half of the time the nuclei 
need to cross each other, matter appears in the CRR. We assume that soon thereafter matter thermalizes and this 
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moment, to, is estimated at about 1 fin feLpl. Then matter starts to expand and cool and after reaching the freeze- 
out temperature it is no longer a thermodynamical system. We denote this moment as t/. .- As we have already 
mentioned in the introduction, this matter is the hadron gas, which consists of all hadrons up to fl~ baryon. The 
expansion proceeds according to the relativistic hydrodynamics equations and for the longitudinal component we have 
the following solution (for details see e.g. ||9|) 

s(t) = — , n B {t) = —j— , (3) 

where sq and n B are initial densities of the entropy and the baryon number respectively. The superimposed 
transverse expansion has the form of the rarefaction wave moving radially inward with a sound velocity c s ■ 
To obtain the time dependence of temperature and baryon number and strangeness chemical potentials one has 

- pd| ) with s, tib and ns given as time dependent quantities. For s(t), ns(t) we 
have expressions (|3|) and ng — since we put the overall strangeness equal to zero during all the evolution (for more 
details see [0). 

The sound velocity squared is given by c\ — ^ and can be evaluated numerically plH . 



to solve numerically equations (2b 



III. J/* ABSORPTION IN HADRONIC MATTER 

In a high energy heavy-ion collision, charmonium states are produced mainly through gluon fusion and it takes 
place during the overlap of colliding nuclei. For the purpose of our model, we shall assume that all cc pairs are created 
at the moment t — 0. Before the fusion, gluons can suffer multiple elastic scattering on nucleons and gain some 
additional transverse momentum in this way fl5|-|l5fl . This manifests for instance in the observed broadening of the 
Pt distribution of J/'J' |T?j] . Following |l7j], we express this effect by the transverse momentum distribution of the 
charmonium states of the form 

5(PT ' eoH ¥i^- exp {"¥i^} ' (4) 

where (Px)"jfq,( e o) is the mean squared transverse momentum of J/^f gained in an A-B collision with the initial 
energy density eo- The momentum can be expressed as (for details see |l8|]) 

(p 2 T)^) = {p 2 T) P j%+K-e, (5) 

with K = 0.27 /m 3 • GeV and {pt) P j% = 1.24GeV 2 taken from a fit to the J/* data of NA38 Collaboration @. 
The expression in (Q) is normalized to unity and is treated as the initial momentum distribution of charmonium states 
here. 

For the simplicity of our model, we shall assume that all charmonium states are completely formed and can be 
absorbed by the constituents of a surrounding medium from the moment of creation. It means that we neglect a 
whole complex process of J/vP formation as presented in |l^,^0|. The main feature of the above-mentioned process is 
that, soon after the moment of production, the cc pair binds a soft gluon and creates a pre-resonance cc — g state, 
from which, after a time of the order of 0.3 fm, a physical charmonium state is formed. This means that the possible 
nuclear absorption of charmonium is, in fact, the absorption of the cc — g state. But the latest has the cross-section 
0ahs — 7.3mfo, which is much higher than J/\f' — Nucleon absorption cross-section ct^jv — 3 — bmb obtained from 
p-A data ||l6| , ^l| , p2"f . This justifies our assumption: taking into account cc — g absorption instead of charmonium 
disintegration in the nuclear matter would only strengthen J suppression. 

According to the above assumption, charmonium states can be absorbed first in the nuclear matter and soon later, 
when the matter appears in the CRR, in the hadron gas. Since these two processes are separated in time, J/^P survival 
factor for a heavy- ion collision with the initial energy density eoi may be written in the form 

Af(eo) =A/" n . m .(e ) -Mh.gXeo) , (6) 

where Af n . m .( e o) an d -^h.g.{ e o) ar e J/^f survival factors in the nuclear matter and the hadron gas, respectively. For 
Nn.m. (eo) we have the usual approximation |^2| Q 

n.m.(eo) = e-np{-a^ N p L} , (7) 
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where po is the nuclear matter density and L the mean path length of the J/^> through the colliding nuclei. For 
the last quantity, we use the expression given in |2^| : 



L(b) = 



1 



%PoTab 



d 2 s T A (s)T b (s-b) T A (s) + T B {s-b) 



(8) 



where Tas(6) = f d 2 s Ta(s)Tb(s — 6), Ta(s) = J dzpA(s,z) is the nuclear density profile function, pa{s,z) the 
nuclear matter density distribution and b the impact parameter. How to obtain eo as a function of b will be presented 
further. 

To estimate Afh.g.( e o) we follow the description presented in ||. We shall focus on the plane z = (z is a collision 
axis) and put J/ty longitudinal momentum equal to zero. Now the pr-dependent J/\l/ survival factor Mh.gXpr) is 
given by (for details see (|]) 



■N'h.gXPT) = I d 2 sf (s,p T )exp 
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where the sum in the power is over all taken species of scatters (hadrons), p v = (E,pr) and q\ = {E' i ^ q) are four 
momenta of J/^ and hadron specie i respectively, v = pr/E is the velocity of the former, <ii states for the absorption 
cross-section of J/^ — hi scattering and v re z,i is the relative velocity of hi hadron with respect to J/'f. When M 
denotes J/^ mass, M = 3097 MeV, iv e ;,i reads 



Vrel.i = ( 1 - 



mfM 2 



(10) 



The upper limit of the time integral in (|^) , tf, is equal to i/. . or to t esc - the moment of leaving by a given J/\l/ 
of the hadron medium, if the final-size effects are considered and t esc < tf .. For <Ji we have assumed that it equals 
zero for {p v + q\) 2 < (2mu + nix) 2 and is constant elsewhere {mo is a charm meson mass, mo = 1867 MeV). For 
hadron specie i we have usual Bose-Einstein or Fermi-Dirac distribution (we neglect any possible spatial dependence 
here) 



mt) = fi(q,t) 



2 Si 



exp 



r El-ma) \ 

\ T(t) / 



(11) 



In the following, we shall consider only J/^P initial distribution /o(s,pr) that factorizes into fo(s)g(pT) and the 
momentum distribution g{pr) will be given by (|J) . We assume at the first step that the transverse size of the 
hadron medium is much greater than tf. a . and also much greater than the size of the area where fo(s) is non-zero. 
Additionally we assume that /o(s) is uniform and normalized to unity. Note that the first assumption overestimates 
the suppression but the second, in the presence of the first, has no any calculable effect here. As a result, Mh.gXPr) 
simplifies to 



Nh.gXPr) = g(pr, eo) • exp • 



dt 



i 

E 
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d 3 q 



7fz(q,t)aiV re i. 



' EE'. 



(12) 



To obtain J\fh. g X e o) one needs only to integrate ( p"2| ) over pr'- 



Nh.gX e o) = / dPT gipr, eo) • exp ■ 



d 3 g 



-fi(q,t)aiV re i. 



Puq t 

* EE' 



(13) 



Now we would like to take the final-size effects and the transverse expansion into account in our model. To do this 
directly, we would have to come back to the formula given by @ and integrate it, instead of ( |l^ ) , over pt- But 
this would involve a five-dimensional integral (the three-dimensional integral over q simplifies to the one-dimensional 
one, in fact) instead of the three-dimensional integral of (|l3|). Therefore, we need to simplify in some way the direct 
method just mentioned above. We shall define an average time of leaving the hadron medium by J/^'s with the 
velocity v produced in an A-B collision at impact parameter b, (t esc )(b,v). Then, we will put this quantity, instead 
of t f. . , as the upper limit of the integral over t in (O) . 
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FIG. 2. View of a Pb-Pb collision at impact parameter b in the transverse plane (z = 0). The region where the nuclei overlap 
has been hatched and its area equals S e //. 

Let us consider an A-B collision at impact parameter b. Since we will compare final results with the latest data of 
NA50 which are for Pb-Pb collisions 0, we focus on the case of A=B here. So, for the collision at impact parameter 
b we have the situation in the plane z = as presented in Fig. |], where S e f / means the area of the overlap of the 
colliding nuclei. We shall assume here, that the hadron medium, which appears in the space between the nuclei after 
they crossed each other also has the shape of S e ff at to in the plane z — 0. And additionally, the transverse expansion 
will start in the form of the rarefaction wave moving inward S e ff at tg- Then, for a J/ty which is at r € S e ff at the 
moment to and has the velocity v we denote by t esc the moment of crossing the border of the hadron gas. It means 
that t esc is a solution of the equation | d + v(t — to) \= Ra — c s (t — to), where Ra = r ■ As (r Q — 1.2 fm) is the 
nucleus radius and d — r — b for the angel between r and v such that the J/^f will cross this part of the edge of the 
area of the hadron gas which was created by the projectile and d = f in the opposite. Having obtain t esc , we average 
it over the angel between r and v, i.e. we integrate t esc over this angel and divide by 27r. Then we average the result 
over S e ff with the weight given by 

T A (r)T B (r-b) 

PjMrl = T AB ib) (W) 

and we obtain (t esc )(b, v). So, the final expression for J\fh. g .( e o) when the transverse expansion is taken into account 
reads 

A^.g.M = J dp T g(pT,eo) ex.pl- J dt^ J J^^fi(^ t ) a i V reLt^r } ■ ( 15 ) 



IV. THE ENERGY DENSITY IN THE CRR 

To compare our theoretical estimations for J/^ survival factor with the experimental data |Q we need to modify 
the latest so as they become eo-dependent instead of E^-dependent (Et is the neutral transverse energy). To do this 
we will use the well-known Bjorken formula 

Ar]b e ffto 

where Arj is the pseudo-rapidity range. Using values of impact parameter b given in we can calculate S e ff for 
each measured Et bin. 

In further considerations we will need the formula for the number of participating nucleons as a function of impact 
parameter b: 

N part (b)= f d 2 s{T A {^)+T B (s-b)} . (17) 
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We found out that the ratio ET/N part is almost constant and for NA50 data varies from 0.28 GeV to 0.31 
GeV with the average value equal roughly to 0.3 GeV and the standard deviation equal to 0.0094 GeV. Therefore we 
assume that for Pb-Pb collisions of NA50 the following approximation is valid: 



Et — 0.3 • N part 



Having put @ into @ we obtain eo as a function of b 



eo(&) = 0.75 



Seff(b) ' 



(18) 



(19) 
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FIG. 3. The initial energy density eo in the CRR for Pb-Pb collisions as the function of impact parameter 6 and for the 
NA50 data of 1996 run |R| 



where we have also used the value A 77 = 1.2 of NA50 |Q. The above function is depicted in Fig. ||. The behaviour 
in low b is the most interesting feature of eo(b). We can see that for b < 7.9 there are two different values b\ and 62 
such that eo(&i) = £0(^2)- Of course, this is the result of direct application of (Cs[) which is some approximation in 
fact. But nevertheless this can suggest that there is a wide range of impact parameter b for which the resulting eo 
is almost constant. Using the dependence between Et and b obtained by NA50 0, we can state the similar for Et, 
i.e. for Et > 40 — 50GeV eo changes weakly with Et- Also the same is true for J/\P suppression what could suggest 
that the quantity of eo reached is the main reason for the suppression. 



V. RESULTS 

To evaluate formulae f^3| ) and ( plf ) we have to know T(t), ^s(i) and jus(i) and how to obtain these functions 
was explained in Sect. U But to follow all that procedure we need initial values so and n B . To estimate initial 
baryon number density n B we can use experimental results for S-S f25| or Au-Au [p6| , p7| collisions. In the first 
approximation we can assume that the baryon multiplicity per unit rapidity in the CRR is proportional to the 
number of participating nucleons. For a sulphur-sulphur collision we have dN B /dy = 6 |25| and 64 participating 
nucleons. For the central collision of lead nuclei we can estimate the number of participating nucleons at 2A — 416, so 
we have dNs/dy = 39. Having taken the initial volume in the CRR equal to ttR\ ■ 1 fm, we arrive at n B = 0.25 fm" 3 . 
This is some underestimation because the S-S collisions were at a beam energy of 200 GeV/nucleon, but Pb-Pb at 
158 GeV/nucleon. From the Au-Au data extrapolation one can estimate n° B = 0.65 fm" 3 |2q] . These values are for 
central collisions, and for the higher impact parameter (a more peripheral collision) the initial baryon number density 
should be much lower. In fact, if we apply the above-mentioned assumption, the initial baryon number density for 
a given collision at the impact parameter b will be proportional to the number of participating nucleons divided by 
S e ff. Therefore, n B (b) will have exactly the same shape as eo(&) presented in Fig. ||. As a result, only for the most 
peripheral collisions n B will be substantially below the value for the central one. So, to simplify numerical calculations 
we will keep n B constant over the all range of b and additionally, to check the possible dependence on n B , we will do 
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our estimations for n g substantially lower, i.e. 



0.05/m 3 . According to our approximation of n B (b), it would 



n B 



and 
and 



correspond to the most peripheral Pb-Pb collisions, b > 14/m. 

Now, to find sq, first we have to solve (|2a| - pc|) with respect to T, us and Mb, where we put e = 
lis — 0. Then, having put T, [is an d [Ib into (pq) we obtain so. Finally, expressing left sides of (2bJ,2cJ) by (| 
after then solving (pb|- pd| ) numerically we can obtain T, /is an d /is as functions of time. In fact, evaluating formulae 
( |l3| ) and ([ll]) we do the following: first, we calculate T = T(t) which turns out to be very well approximated by the 
expression 



T(t) T • r a 



(20) 



and then we put this approximation into (13) and ( |l5| ) . And for fJ,s(t) an d (J-B(t) in fj (q. t) we put solutions of 
( ph| , |2c| ) where given by (||) and ns = and T is given by (^0|). But the exponent a in ( f20| ) has proven not to be 
unique for the whole range of To considered here. One gets different values of the initial energy density eo for different 
values of the impact parameter b and for different geometry of the collision process. So b dependent a gives also b 
dependent freeze-out time tf. . The formula (|l9| ) is calculated for different values of b and Eqs. (||a| - |2c|) are solved. 
We have evaluated the suppression factor up to eo = 3.5GeV"//ra 3 . This gives e.g. the maximal possible To, To iTnax , 
equal to 219.3MeV ( for n% = 0.65 fm~ 3 ), 22AMeV (for n% = 0.25/m~ 3 ) or 224.8MeV (for n° B = 0.05/m- 3 ). ' 

This procedure allows to evaluate J/ iff survival factor given by (O) . Because of the lack of data, we shall assume 
only two types of the cross-section, the first, Ub, for J/vf-baryon scattering and the second, a m , for J/vP-meson 
scattering. For o\> we put <7& = (Jj/^n- As far as a m is concerned, we assume that this cross-section is 2/3 of the 
corresponding cross section for baryons, which is due to the quark counting. In the following, we will use values of 



J /\P — Nucleoli absorption cross-section a 



J/ipN 



3 — hmb obtained from p-A data p 



At the beginning, to 



illustrate how the value of power a influences J/^> suppression we present in Fig. |l| two results: the first for a = | 
(which is the exact value for a free massless gas) and the second for a = ^ (which is the approximate value for the 
hadron gas and Tq = 200MeT^). We can see that the suppression improves more than twice for the highest eo indeed. 
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FIG. 5. Same as Fig. | but for7>. . 
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FIG. 4. J/* suppression in the longitudinally ex- FIG. :,. San,,. I-V I l>nt forT, „ - M) \Jr\\ 

panding hadron gas with the " infinite" transverse size 
and for n% = 0.25 /m -3 and T f .„. = 140 MeV: a) 
(Tb — 3 mb, o m — 2 mb; b) Ob = 4 mb, a m = 2.66 mb; 
c) <Tb — 5 mb, o m — 3.33 mb. Triangles and crosses 
represent the S-U data of NA38 and the Pb-Pb data 
of NA50 Collaborations respectively. 

To make our investigations more realistic we have to take into account that only about 60% of J/^S measured are 
directly produced during collision. The rest is the result of \ (~ 30%) and ip' (~ 10%) decay [J2y] . Therefore the 
realistic J/^> survival factor should read 



•AA(eo) = 0.6A/>/^(eo) + 0.3A^(e ) + 0.17V^(e ) 



(21) 
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where A/j/v>( e o), N x (eo) and -^/j'( £ o) are given also by formulae ([| - |l5| ) but with 

cj/rpN — (J j/ipN,OxNi (T il>'N and A/ = 



Kj/i>,K x , K,i 



respectively. The remaining problem is whether formula 
and they shows that (jo>r)™ Pb 



Mj^,M x ,M^ 



is valid for x and tp'. There are data for (jpj)% 



So, we assume that the above is also true for (£4)^/ B (e), i.e. 



with (p 2 T )j^(e) given by (|) 



(p|)^(e)=1.4{p|)^(e) 
For x we believe that the inequality 



(P 2 T ) 



2 \AB 



,//* 



< 



(Pt) 



2 \AB 



< 



2 \ AB 



(22) 



(23) 



should be valid and therefore assume that (E 



is true also in this case. Anyway, the exact form of {Pr) x B {e) or 
(Pt)^ /B ( £ ) i s n °t very important because we checked that the suppression depends on this form very weakly. First, 
we put Kji^) = and the resulting J survival factor (for direct J/'J's) differs only a few percent for the highest 
eo from the one calculated with formula (||) unchanged. Second, when we use expression (||) also for x and ip', the 
evaluated suppression factor is the same as that calculated with the use of (E2), as far as plots are concerned. 
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FIG. 6. J/\& suppression in the longitudinally 
expanding hadron gas with the "infinite" trans- 
verse size and for at — 4 mb, a m = 2.66 mb 
and Tf.o. = 140 MeV. The curves correspond 



to n B = 0.05 /m~ J (solid), 



= 0.25 fm' 
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FIG. 7. Same as Fig. | but for !>.„. = 100 MeF. 



(short-dashed) and n° B = 0.65 fm~ 3 (dashed). Tri- 
angles and crosses represent the S-U data of NA38 
and the Pb-Pb data of NA50 Collaborations Q , re- 
spectively. 

To complete our estimations we need also values of cross-sections for x ~ baryon and ip' — baryon scatterings (we 
will still hold that xWO — rneson cross-section is | of x(V'') — baryon cross-section). Since J/^f is smaller than x or 
-0', x ~ baryon and ip' — baryon cross-sections should be greater than J — baryon one. For simplicity, we assume 
that all these cross-sections are equal. This means that we underestimate J /\& suppression, here. The final results of 
calculations of ( |l3"l) are presented in Figs. ^0for various sets of parameters of our model (which are T/. ., n^, Cb). We 
performed these calculations for two values of T/. c . = 100, 140MeV which agree fairly well with values deduced from 
hadron yields |pq| . For comparison, also the experimental data are shown in Figs. The experimental survival 
factor is defined as 



(24) 



8 



where 



, AB(pp) 

AB(pp) 
DY 



is the ratio of the J/ty to the Drell-Yan production cross-section in A-B(p-p) interactions times 



the branching ratio of the J/*? into a muon pair. The values of the ratio for p-p, S-U and Pb-Pb are taken from [Q. 
Note that since the equality a^y = oV dy ' has been confirmed experimentally up to now |^8| , formula ( |2~4| ) reduces 
to 



AC 



ABo% Y 



which is also given as the experimental survival factor, for instance, in 



(25) 
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FIG. 8. J/\P suppression in the longitudinally and 
transversely expanding hadron gas for the uniform 
nuclear matter density distribution and aj, = 4 mb, 
cr m = 2.66 mb and T>.„. = 140 MeV. The curves 
correspond to n° B = 0.25 fm~ 3 , c B — 0.45 (solid) and 
n B = 0.65 fm~ 3 , c s = 0.46 (dashed). Triangles and 
crosses represent the S-U data of N A38 and the Pb-Pb 
data of NA50 Collaborations |Q respectively, but the 
data are "cleaned out" from the contribution of J/$ 
scattering in the nuclear matter, in accordance with 
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FIG. 9. J/ty suppression in the longitudinally 
and transversely expanding hadron gas for the 
Woods-Saxon nuclear matter density distribution and 
a b = 4 mb, a m = 2.66 mb and T/. D . = 140 MeV. The 
curves correspond to n% = 0.25 fm~ 3 , c s = 0.45 



(solid) and n° B — 0.65 fm 



0.46 (dashed). 



Triangles and crosses represent the S-U data of NA38 
and the Pb-Pb data of NA50 Collaborations Q re- 
spectively, but the data are "cleaned out" from the 
contribution of J/*? scattering in the nuclear matter, 
in accordance with ( pi] ). 



Coming back to examination of our first results presented in Figs. we can see that the most of Pb-Pb data 
are below the region of suppression obtained for chosen values of parameters in our model. Generally, the theoretical 
curves decrease slower with increasing eo, whereas the data show rather abrupt fall just above eo = 2GeV/ fm 3 . Note 
that the dependence on the initial baryon number density is substantial but for higher values of n° B , rather. The 
lower the initial baryon number density, the deeper the suppression. There are two reasons for such a behaviour: the 
first, for the higher baryon number density, there are less non-strange heavier mesons p, lu in the hadron gas of the 
same eo, but these particles create the most weighty fraction of scatters, for which reaction ([!]) have no threshold 
at all; the second, the freeze-out time tf. . decreases with increasing mP B for a given eo in our model. For instance, 
for e = 3.5GeV/fm 3 and T/. , = 140MeP we have a = 0.172, 0.175, 0.183 and f/. . = 15.7, 14.7, 11.6/m for 
n B = 0.05, 0.25, 0.65/m~ 3 , respectively. We can see also that the value <Jb = 3mb is too small to obtain results 
comparable with the data, so we will leave aside this value in further investigations. 

Now we will include the finite-size effects into our model, i.e. we will take into account that the realistic hadron 
gas has a finite transverse size. This will be done in form of the rarefaction wave moving inward S e ff with the 
sound velocity c s . How to obtain this velocity has been mentioned in Sec. 2 (see als o [12| ). With the finite-size effects 
included, the final expression for J/\l/ survival factor Nh.g.( e o) w iU be given by (15) . To make our investigations 
much more realistic we will also include the possible J/^ disintegration in nuclear matter, which should increase J/^ 
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suppression by about 10% pi"| . But to draw also S-U data in figures, instead of multiplying Nh.g. by M n . m . given by 
([7|), we divide Af eX p by appropriate Af n . m ., i-e. we define "the experimental J/ty hadron gas survival factor" as 



1 .3 "T 

1 .2 - 




0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 - 

0.0 1 1 1 1 1 1 1 

12 3 

£„ [GeV/fm 3 ] 

FIG. 10. Same as Fig. ^| but for <jj, = 5 mb, a m = 3.33 mb. 

Mex-p = CXp {(Tj/tPNPqL} ■ Mexp ■ (26) 

and values of this factor are drawn in Figs. |§|-|lC| as the experimental data. 
At the beginning we will consider a uniform nuclear matter density 

p A (s, z) = PA(r) = { A. = (f '-o 3 ) 1 , \f \ZRA . (27) 
10 , \r\> R A 

The results of numerical estimations of (|l5|) and ( p6| ) are depicted in Fig. || for two values of the initial baryon 
number density, n° B = 0.25, 0.65 fm~ 3 . The curve for n° B = 0.05 fm~ 3 almost covers the curve for n° B = 0.25 fm~ 3 
and the maximal difference between these curves do not exceed 1.9% (which is the difference for the highest possible 
eo, i-e. e ~ 2GeV/ fm~ 3 ), so for clearness of the figure we do not draw it. The two values of the speed of sound 
are the maximal values of this quantity possible in the range [T/. D . = 140 MeV, To tma x] for the above-mentioned two 
cases of n B . In fact, we have checked that the results almost do not depend on c s (allowed in the range) and the 
difference (seen only for the quantity eo) between J/^> survival factors for the maximal and the minimal values of 
c s in the range [TV . = 140 MeV,To, max ] are less than 0.6%. Note that the theoretical curves in Figs. || (the same 
will happen in Figs. |9|-|To|) are two-valued around eo = 2GeV/ fm~ 3 . This is the result of our approximation of e (&) 
given by (|l9|). This expression allows for two different values of 6, which give the same eo in some range of the impact 
parameter (see Sec. 4). 

It has turned out also that in the case of the transverse expansion, the results almost do not depend on the T/ . (for 
Xy o . £ [100, 140}MeV). And the maximal difference (seen for the quantity eo) between curves for T/. . = lOOMeV 
and Tf. . = UOMeV do not exceed 2.7% (n% = 0.05/m" 3 ), 2.6% (n% = 0.25/m~ 3 ). This is because the freeze-out 
time resulting from the transverse expansion, t/. .,trona — Ra/c s (if we assume a central collision and c s constant), is 
of the order of the freeze-out time resulting from the longitudinal expansion for T/. D . = 140M eV . Namely, for Pb and 
c s = 0.45 we have i/. .,trans — 15.8/m which is very similar to values of tf. a , for T/ . = 140Mey given earlier. For 
Tf. . = WQMeV, t f . . = 111.0, 101.0, 72.5/m for n° B = 0.05, 0.25, 0.65/m- 3 respectively, so the hadron gas ceases 
because of the transverse expansion much earlier. 

We repeated our estimations of formula ( [j"5"| ) also for the more realistic nuclear matter density distribution, namely 
for the Woods-Saxon distribution with parameters taken from [[30|. The results are presented in Figs. p[jltj|. We can 
see that these curves fit the data a little bit better than those obtained within the assumption of the uniform nuclear 
matter density distribution (cf. Fig. || and Generally, taking into account also the transverse expansion changes 
the final (theoretical) pattern of J/^ suppression qualitatively. First of all, the curves for the case including the 
transverse expansion are concave (what is more clearly seen for lower n B ) as the data suggest should be, in opposite 
to the case with the longitudinal expansion only, where curves are convex. But still, theoretical curves are not steep 
enough to cover the data area completely. We would like to note at this point that there is some ambiguity in 
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calculation of (Et) and b in NA50 experiment, since the range of eo obtained from ( Jlq ) with the use of (Et) and b 
for 1995 and 1996 Pb-Pb runs are different. For 1995 run fjlj we have e S 2.6 - 2.9GeV/fm- 3 and for 1996 run § 
eo — 2.0 — 2.1GeV/ fm~ 3 , both estimates are for (Et) > 40GeV. Of course, where exactly the data points should be 
placed is crucial for the valuation of the correctness of the shape of theoretical curves. 

When the additional disintegration in the nuclear matter is included, also the magnitude of J suppression is 
comparable with the data, but rather for greater ov But note that since we have one overall charmonium-baryon 
cross-section ab, our final results underestimate the suppression (for x~ > "0' — baryon scattering the cross-section 
should be greater than for J/^f). 

As a final remark, we think that it is difficult to exclude J/^f scattering in the hot hadron gas entirely, as the 
reason for the observed J /\& suppression at this point. In our model the most crucial parameter is the charmonium- 
baryon inelastic cross-section and the final results depend on its value substantially. Therefore it is of the greatest 
importance to establish how this cross-section behaves in the hot hadron environment. Some work has been done into 
this direction p2|-p4|, but results presented there differ from each other and are based on different models. 
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